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Abstract. We discuss the motivation for resummation of the effects of initial- 
state soft gluon radiation, to all orders in the strong coupling strength, for pro¬ 
cesses in which the near-threshold region in the partonic subenergy is important. 
We summarize our calculation of the total cross section for top quark produc¬ 
tion at hadron colliders. Comments are included on the differences between our 
treatment of subleading logarithmic terms and other methods. 


INTRODUCTION AND MOTIVATION 


In inclusive hadron interactions at collider energies, tt pair production pro¬ 
ceeds through partonic hard-scattering processes involving initial-state light 
quarks q and gluons g. In lowest-order perturbative quantum chromodynam¬ 
ics (QCD), at 0(al), the two partonic subprocesses are q + q —> t + t and 
g + g —> t +1. Calculations of the cross section through next-to-leading order, 
0(cXs), involve gluonic radiative corrections to these lowest-order subprocesses 
as well as contributions from the q + g initial state [1]. In this paper, we de¬ 
scribe calculations that go beyond fixed-order perturbation theory through 
resummation of the effects of gluon radiation [2-4] to all orders in the strong 
coupling strength a s . 

The physical cross section is obtained through the factorization theorem 


aij(S,7n ) 


4 m 2 

~S~ 



dV^ijirh ri&ijirhm, (i). 


( 1 ) 


The square of the total hadronic center-of-mass energy is S, the square of the 
partonic center-of-mass energy is s, m denotes the top mass, // is the usual 
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factorization and renormalization scale, and $^(77, 77) is the parton flux. The 
variable 77 = — 1 measures the distance from the partonic threshold. The 

indices ij G {qq, gg} denote the initial parton channel. The partonic cross 
section by,-(77, 777,77) is obtained either from fixed-order QCD calculations [ 1 ], 
or, as described here, from calculations that include of resummation [ 2 - 4 ] to 
all orders in a s . We use the notation a = a(g = m ) = a s (m)/ 7 i. The total 
physical cross section is obtained after incoherent addition of the contributions 
from the the qq and gg production channels. 

Comparison of the partonic cross section at next-to-leading order with its 
lowest-order value reveals that the ratio becomes very large in the near¬ 
threshold region. Indeed, as 77 —> 0 , the “K- -factor” at the partonic level K{rf) 
grows in proportion to a In 2 (77). The very large mass of the top quark notwith¬ 
standing, the large ratio K(rj) makes it evident that the next-to-leading order 
result does not necessarily provide a reliable quantitative prediction of the 
top quark production cross section at the energy of the Tevatron collider. 
Analogous examples include the production of hadronic jets that carry large 
values of transverse momentum and the production of pairs of supersymmetric 
particles with large mass. 


GLUON RADIATION AND RESUMMATION 

The origin of the large threshold enhancement may be traced to initial- 
state gluonic radiative corrections to the lowest-order channels. We remark 
that we are calculating the inclusive total cross section for the production of 
a top quark-antiquark pair, i.e., the total cross section for t + t + anything. 
The partonic subenergy threshold in question is the threshold for t + t+ any 
number of gluons. This coincides with the threshold in the invariant mass of 
the t + t system for the lowest order subprocesses only. 

For i + j —> t + t + g, we define the variable z through the invariant 
(1 — z) = where k and pt are the four-vector momenta of the gluon and 
top quark. In the limit that z —> 1 , the radiated gluon carries zero momentum. 
After cancellation of soft singularities and factorization of collinear singulari¬ 
ties in 0 (ag), there is a left-over integrable large logarithmic contribution to 
the partonic cross section associated with initial-state gluon radiation. This 
contribution is often expressed in terms of “plus” distributions. In it is 

proportional to a 3 ln 2 (l — z). When integrated over the near-threshold region 
1 > z > 0, it provides an excellent approximation to the full next-to-leading 
order physical cross section as a function of the top mass. The goal of gluon 
resummation is to sum the series in a n+2 ln 2n (l — z) to all orders in a in order 
to obtain a more trustworthy prediction. 

Different methods of resummation differ in theoretically and phenomeno¬ 
logically important respects. Formally, if not explicitly in some approaches, 
an integral over the radiated gluon momentum 2 must be done over regions 



in which z —» 1. Therefore, one significant distinction among methods has to 
do with how the inevitable “non-perturbative” region is handled. 

The method of resummation we employ [3] is based on a perturbative trun¬ 
cation of principal-value (PV) resummation [5]. This approach has an im¬ 
portant technical advantage in that it does not depend on arbitrary infrared 
cutoffs. Because extra scales are absent, the method permits an evaluation of 
its perturbative regime of applicability, i.e., the region of the gluon radiation 
phase space where perturbation theory should be valid. We work in the MS 
factorization scheme. 

Factorization and evolution lead directly to exponentiation of the set of 
large threshold logarithms in moment (n) space in terms of an exponent E px . 
The function E px is finite, and lim^-^ E PV (n, m 2 ) = — oc. Therefore, the 
corresponding partonic cross section is finite as z —> 1 (n —> +oo). The 
function E px includes both perturbative and non-perturbative content. The 
non-perturbative content is not a prediction of perturbative QCD. We choose 
to use the exponent only in the interval in moment space in which the pertur¬ 
bative content dominates. We derive a perturbative asymptotic representation 
of E(x,a(m )) that is valid in the moment-space interval 

1 < x = Inn < t = —- . (2) 

2 ab 2 y ’ 

The coefficient b- 2 = (11 Ca — 2n/)/12; the number of flavors nj = 5; C qq = 
Cf = 4/3; and C gg = Ca = 3. 

The perturbative asymptotic representation is 

N(t)+1 p +1 

Eij(x, a) ~ Eij(x,a, N(t)) = 2C i:j ^ a p ^ s^ p x 3 . (3) 

P =1 3 =0 

Here 

= -f-r't-l)' ,+J 2% +1 -,(p - 1 )!/j! ; (4) 

and T(1 + z) = J2kLo c k zk , where T is the Euler gamma function. The num¬ 
ber of perturbative terms N(t) in Eq. (3) is obtained [3] by optimizing the 

asymptotic approximation E(x,a) — E(x,a, N(t)) — minimum. Optimiza¬ 
tion works perfectly, with N(t) = 6 at m = 175 GeV. As long as n is in 
the interval of Eq. (2), all the members of the family in n are optimized at 
the same N(t), showing that the optimum number of perturbative terms is a 
function of t, i.e., of m only. 

Resummation is completed in a finite number of steps. Upon using the 
running of the coupling strength a up to two loops only, all monomials of the 
form a k ln fc+1 n, a k ln fc n are produced in the exponent of Eq. (3). We discard 
monomials a k ln fc n in the exponent because of the restricted leading-logarithm 



universality between tt production and massive lepton-pair production, the 
Drell-Yan process. 

The exponent we use is the truncation 

N(t )+1 

E i:j (x,a,N) = 2C i:j ® P SpX p+l , (5) 

p =i 

with the coefficients s p = s p+ i tP = & 2 _ 1 2 p /p(p + !)• This expression contains 
no factorially-growing (renormalon) terms. It is valuable to stress that we 
can derive the perturbative expressions, Eqs. (2), (3), and (4), without the 
principal-value prescription, although with less certitude [3]. 

After inversion of the Mellin transform from moment space to the physically 
relevant momentum space, the resummed partonic cross sections, including all 
large threshold corrections, can be written 

a% pert (i h m) = f ZmaX m, z). ( 6 ) 

“ Zmin 

The leading large threshold corrections are contained in the exponent 
Eij(x,a), a calculable polynomial in x. The derivative < 7 ^( 77 , m, z) = 

d(d’ij\r)i r n,z))/dz, and df/ is the lowest-order O(a^) partonic cross section 
expressed in terms of inelastic kinematic variables. The lower limit of inte¬ 
gration, z m i n , is fixed by kinematics. The upper limit, z rnax < 1, is set by the 
boundary between the perturbative and non-perturbative regimes, well spec¬ 
ified within the context of our calculation. Its presence assures us that our 
integration over the soft-gluon momentum is carried out only over a range in 
which poorly specified non-universal subleading terms would not contribute 
significantly even if retained. We cannot justify continuing the results of 
leading-logarithm resummation into the region z > z max . 

Perturbative resummation probes the threshold down to 77 > r] 0 — (1 — 
Zmax)/ 2- Below this value, perturbation theory is not to be trusted. For m 
= 175 GeV, we determine that the perturbative regime is restricted to values 
of the subenergy greater than 1.22 GeV above the threshold (2m) in the 
qq channel and 8.64 GeV above threshold in the gg channel. The difference 
reflects the larger color factor in the gg case. The value 1.22 GeV is comparable 
to the decay width of the top quark, a natural definition of the perturbative 
boundary and by no means unphysically large. 

PHYSICAL CROSS SECTION 

Other than the top mass, the only undetermined scales are the QCD fac¬ 
torization and renormalization scales. We adopt a common value p for both. 
In Fig. 1, we show our total cross section for ff-production as a function of 
top mass in pp collisions at a /S = 1.8 TeV. The central value is obtained 



with the choice /x/m = 1, and the lower and upper limits are the maximum 
and minimum of the cross section in the range /i/m e {0.5,2}. At m = 175 
GeV, the full width of the uncertainty band is about 10% . As is to be ex¬ 
pected, less variation with /x is evident in the resummed cross section than 
in the next-to-leading order cross section. In estimating uncertainties, we do 
not consider explicit variations of our non-perturbative boundary, expressed 
through z max . This is justified because, for a fixed m and p, z max is obtained 
by enforcing dominance of the universal leading logarithmic terms over the 
subleading ones. Therefore, z max is derived and is not a source of uncertainty. 
At fixed m, the boundary necessarily varies as p and thus a vary. 

Our calculation is in agreement with the data [6]. We find a tt (rn = 
175 GeV, y/S = 1.8 TeV) = 5.52/^42 pb. This cross section is larger than 
the next-to-leading order value by about 9%. The top quark cross sec¬ 
tion increases quickly with the energy of the pp collider. We determine 
< 7 a (m = 175 GeV, y/S = 2 TeV) = 7.56io;55 pb. The central value rises 
to 22.4 pb at y/S = 3 TeV and 46 pb at y/S = 4 TeV. 

Extending our calculation to much larger values of m at = 1.8 TeV, we 
find that resummation in the principal qq channel produces enhancements over 
the next-to-leading order cross section of 21%, 26%, and 34%, respectively, for 
m = 500, 600, and 700 GeV. The reason for the increase of the enhancements 
with mass at fixed energy is that the threshold region becomes increasingly 
dominant. Since the qq channel also dominates in the production of hadronic 
jets at very large values of transverse momenta, we suggest that on the order 
of 20% of the excess cross section reported by the CDF collaboration [7] may 



FIGURE 1. Inclusive total cross section for top quark production. The dashed curves 
show the upper and lower limits while the solid curve is our central prediction. CDF and 
DO data are shown. 




be accounted for by resummation. 


OTHER METHODS OF RESUMMATION 


Two other groups have published calculations of the total cross section 
at m = 175 GeV and y/s = 1.8 TeV: cr“(LSvN [2]) = 4.951 q;4o pb; and 
<r tt (CMNT [4]) = 4.751^68 pb- From a numerical point of view, all agree 
within their estimates of theoretical uncertainty. However, the resummation 
methods differ as do the methods for estimating uncertainties. Both the cen¬ 
tral value and the band of uncertainty of the LSvN predictions are sensitive to 
their infrared cutoffs. To estimate theoretical uncertainty, we use the standard 
/i- variation, whereas LSvN obtain theirs primarily from variations of their cut¬ 
offs. It is difficult to be certain of the central value and to evaluate theoretical 
uncertainties in a method that requires an undetermined infrared cutoff. 

The group of Catani, Mangano, Nason, and Trentadue (CMNT) [4] calcu¬ 
late a central value of the resummed cross section (also with g/m = 1) that 
is less than 1% above the exact next-to-leading order value. There are simi¬ 
larities and differences between our approach and the method of CMNT. We 
use the same universal leading-logarithm expression in moment space, but dif¬ 
ferences occur after the transformation to momentum space. The differences 
can be stated more explicitly if we examine the perturbative expansion of the 
resummed hard kernel 7ifj(z,a). If, instead of restricting the resummation 
to the universal leading logarithms only, we were to use the full content of 
we would arrive at an analytic expression that is equivalent to the 
numerical inversion of Ref. [4], 


Tfg ~ 1 + 2a C, 


y 


In 2 (1 - z) + 2y E ln(l - z) 


+ 0(a 2 


(7) 


In terms of this expansion, in our work we retain only the leading term ln 2 (l — 
z) at order a, but both this term and the non-universal subleading term 
27 #ln(l — z) are retained in Ref. [4]. If this subleading term is discarded 
in Eq. (7), the residuals Sij/a^ LO defined in Ref. [4] increase from 0.18% to 
1.3% in the qq production channel and from 5.4% to 20.2% in the gg channel. 
After addition of the two channels, the total residual S/a NLO grows from the 
negligible value of about 0.8% to the value 3.5%. While still smaller than 
the increase of 9% that we obtain, the increase of 3.5% vs. 0.8% shows the 
substantial influence of the subleading logarithmic terms retained in Ref. [4], 
We judge that it is preferable to integrate over only the region of phase 
space in which the subleading term is suppressed numerically. Our reasons 
include the fact that the subleading term is not universal, is not the same 
as the subleading term in the exact 0(a 3 ) calculation and can be changed if 
one elects to keep non-leading terms in moment space. The subleading term 
is negative and numerically very significant when it is integrated throughout 



phase space (i.e., into the region of z above our z max ). In the qq channel at 
m = 175 GeV and y/S = 1.8 TeV, its inclusion eliminates more than half of 
the contribution from the leading term. Although the goal is to resum the 
threshold contributions responsible for the enhancement of the cross section at 
next-to-leading order, the method of Ref. [4] does not reproduce most of this 
enhancement. The influence of non-universal sub leading terms is amplified at 
higher orders where additional subleading structures occur in the approach of 
Ref. [4] with significant numerical coefficients proportional to n 2 , £(3), and 
so forth. In our view, the results of a leading-logarithm resummation should 
not rely on subleading structures in any significant manner. The essence of 
our determination of the perturbative boundary z max is precisely that below 
z max subleading structures are also numerically subleading, whether or not 
classes of poorly substantiated subleading logarithms are included. For a more 
detailed discussion of these and other points of difference, including issues of 
momentum conservation and factorial growth, see section VI of Ref. [8]. 

Since the large threshold logarithms are mastered by resummation, the the¬ 
oretical reliability of the resummed result is greater than that of a hxed-oder 
calculation. Our analysis and the stability of our cross section under varia¬ 
tion of the hard scale fi provide confidence that our perturbative resummation 
yields an accurate calculation of the inclusive top quark cross section at Teva- 
tron energies. 
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